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ABSTRACT 


Digital filter design has become instrumental in many 
fields of electronics. The objectives of this work 1s to 
Study the benefit of using computer design tools to find 
equivalent digital filter realizations of analog filters. 
An analog video bandpass filter for TV Intermediate 
Frequency signal extraction was selected as an example. 

First we introduce the Finite Impulse Response Filter 
(FIR) and the Remez program that is used for optimal FIR 
filter design. Other computer design tools such as ILS are 
used to design a number of existing bandpass filters. 
Emphasis was placed on the use of zero-one decomposition of 
filters to eliminate the number of multipliers and thus 


reduce hardware requirements. 
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TNT RODUCTION 


Digital filter design is important in many electronic 
fields. With advances in digital circuits, digital filter 
design can be applied to many new areas. What kind of tools 
are the best and what kind of methodology can be used to 
design filters effectively are the issues addressed in this 
thesis. First, the Finite Impulse Response Filter (FIR) and 
the Remez program for optimal FIR filter design are 
introduced. Then, existing computer aided design tools are 
used to design several bandpass filter examples. A 
multi-band analog video bandpass filter for TV Intermediate 
Frequency (IF) signal extraction was selected as an example. 
The objective of this work is to study the benefit of using 
Cemeweecmacsion tools to find an equivalent digital filter 
that 1S Similar to the analog one. In terms of hardware, it 
1S 1lmportant to examine the implementation issues of the 
digital filter as well, with interest concentrated on 
zero-one decomposition of a filter to eliminate using 


Mmaleti piers. 


ne BACKGROUND OF FIR FILTER AND CHARACTERISTICS 

The design of nonrecursive FIR digital filters is based 
eae enecmerollowing relationships “"PRe®. 1]. For a causal 
nonrecurSive system the FIR filter can be described by the 


Girference equation, 


k=L 

a ee eas alenica) 
K=O 

(1.1) 

jné 


lt is. found sella tant we ehem npUteT seine then the steady- 


state system output is, 


y (n) = eJn® Hy (e)9) (iva 


where, H(e I8) is the system frequency response. EXDanaame 
; , jno 
the right-side of equation (1.1) for the same input <u ee 


yields the following steady-state output 


née 18 in8 -+ 
e7 H (2? ) = be jn as aoe ae ere 
0 1 
coe et bee ase 
(i230 
— , jnt 
Eliminating e yields, 
, 9 n=L _- | 
H fae ) = h (n) e jne 
n=0 
(1.4) 


where h(n)=by is the impulse response. 

For a casual filter, h(nm)=0 for n<0, with a fine 
integer of duration L. Equation (1.4) establishes a direct 
relationship between the impulse response of a FIR filter 
and the system transfer function. 

The goal of any design is to determine the filter 


Goefficients (Or wetoqhesi = ba, San, eee »Dyn such that a 


16 
desired frequency response characterl1stic H(e? ie cane Oe 
achieved. But, why 1s an FIR digital filter chosen instead 
of other forms of digital filter? The advantages of the FIR 
filter can be described in the following. [Ref. 2] 

1. FIR filters with exact linear phase can be easily 
designed. Linear phase filters are important for 
applications where frequency dispersion due to a 
fommuemmearl phase Shitt 2s harmful, e.g., video and 
speech signal processing. 


Zee ecasrency  Feakryzea monrecursively by direct 
convolution, are always stable. 


The possible disadvantages of FIR filters are: 

1. A large value of N, the impulse response duration, is 
required to approximate sharp cutoff filters. Hence, 
an expensive implementation is required to realize 
Ssulenm- L£llters. 

2. The delay of linear phase FIR filters may not be an 
integer number of samples. This noninteger delay can 
lead to problems in signal processing applications. 
Consequently, a filter with an odd number length to 
yield integer delay is usually preferred. 

Be LAP Seen DALIONS OF FIR FILTER 

There are several methods for hardware realization of an 
FIR filter. The goal is to speed up the execution and make 
possible "real time" simulations in desired digital systems 
with lowest costs. In this section various digital hardware 
realizations of digital filters including direct and cascade 
FIR realizations are discussed. 

1. Direct Form FIR Hardware 

Smemruealizagcelonmeror a direct form FIR filter is 

shown in Figure 1. Figure 2 shows a simple structure for 


realizing the filter using a single computational element 


3 


Ret «.V2in It consists Of a MUMeIplicuwand sonmeceo ms A 
Shift register is used to hold the filter states, andvamee 
is used for the coefficients. By means of a multiplier and 
an accumulator, a single output sample can be computed by 


Successive addition as the shift register content circu- 
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Figure 1 Direct form FIR filter. [Ref. 2] 
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Figure 2 Structure for implementation of direct 
form FIR filter. [Ref. 2] 
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2. Cascade Form FIR Hardware 
Figure 3 shows an example of a three section filter, 
Sxemescemeren Delng FOUREh order. Figure 1 shows the direct 
implementation of (1 + a,271 + aoZ74 +....4+ ay 92712). 
Figure 3 is the cascade implementation of the factored form 
of the above the equation, 1l.e., 
Gee ees baz ty (i = cy27t +....+ ©4274) 
(ee ai z-+ ee. + daz-4) 
The implementation in Figure 3 needs multipliers and adders. 
Sometimes, it is difficult to achieve real time response in 
this implementation due to the speed limitation of the 
mreebp lier . 
How to eliminate the need for multiplication in a 
bandpass filter design? It is possible to accomplish this 


objective using zero-one decomposition methodology [Ref. 3]. 


os 
> 
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Figure 3 Cascade FIR filter. 


Figure 4 shows an implementation of a bandpass 
Pr ceery 


1 - 227-1 - 2-2 + 527-3 - 227-4 - 52-5 + 427-5 + 4277 - 
ieee ae re Oe ell S b7-12 4 7-13 (1.5) 


> 


This bandpass filter has frequency response similar to the 
desired transfer function for color signal used ina 
receiver system. The same transfer function can be written 
as: 


(l= 272) Sy ee ee ae 


( 156) 
Equation (1.6) shows a special way to decompose the transfer 
function where all coefficients are +1,-1 or 0. If the 
desired filter transfer function can be decomposed into this 
kind of form, both processing time and hardware implementa- 
tion cost can be saved. For example, according to the 
expression in (1.5), implementing it in direct form required 
12 multipliers and 13 adders. Implementing the decomposed 
form in (1.6) required only 6 adders. What procedure can be 
used to obtain this decomposed expression? Is this the best 
way to simulate an original filter? These are the issues 
discussed in the Chapter IV. In Chapter I1, optimizatvon 
principle and the characteristics of the Remez program are 
Summarized. In Chapter III cascade and parallel implementa- 
tion of the design example mentioned in Chapter II will be 
discussed. Finally all studies are summarized in Chapter 


IV. 
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decomposition that eliminates the multiplication inherent 


in the filter implementation. [Ref. 3] 
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In this chapter the characteristics of the computer 
tools for designing FIR digital filters are discussed. The 
uses of this computer tool and the advantages of using this 
tool for digital filter design are the important issues 
addressed here. Using the video Intermediate Frequency (IF) 
analog filter as an example we investigate the necessary 


steps to implement an equivalent digital filter. 


A. COMPUTER-AIDED DESIGN USING THE REMEZ EXCHANGE ALGORITHM 
Optimization techniques are used widely in operation 
research, economics, and other related fields. However, 
only recently have such techniques gained acceptance as 
tools for the design and evaluation of electronic circuits. 
1. Optimal Filter Principle 
Since we are only concerned with FIR filter design, 
the Remez exchange algorithm will be chosen as the computer 
aided design tool. The following discussion consists of 
excerpts from the book of Theory and Application of Digital 
Signal Processing [Ref. 2] and the thesis by J. V. England. 
(Reta i). The basis of this computer aided design (CAD) 
technique is optimization. A desired filter frequency 
response is approximated by a particular filter whose 


coefficients are to be determined. The accuracy of the 


adoeeseexttcaeronm 1s evaluated according tO some criterion, 
usually an error function, that indicates how large a 
disparity exists between the desired filter frequency 
response and the approximating filter frequency response. 
Variable parameters of the approximating function are then 
"adjusted" to optimize the filter design in terms of this 
Crisverion. In other words, the solution is optimal (in the 
sense that the peak approximation error over the entire 
interval of approximation 1S minimized) and unique. And the 
mathematical basis for this algorithm is the weighted 
Chebyshev approximation. 

A summary of the approximation and error functions 
port ecmamgorithmetol lows. It has been shown in (2.1) 
(Ref. 2] that the frequency response for the four cases of 
linear phase filter, i1.e., even or odd symmetry with an 
even or odd number of terms, can be written in the form: 

H (el) = ea (nl) / 263) (7 / 2) H (el) 8) (2h 
where ' (eJ®) is a real-valued function used to approximate 
the desired filter's magnitude specifications and the 
remaining terms approximate the desired phase. 

A measure of how well the designed filter frequency 
response approximates the desired filter frequency response 
is mequamed. The weighted Chebyshev approximation uses an 


error function defined as follows: 


E (6) = W (8) [Hp(el®) - H (e5%)] (2.2) 


where, 


Hp (e]®) = the desired frequency response 
H (eJ°) = the designed frequency response 
W (8) = weighting function 


Eel = lennon function 

Thus, the Chebyshev approximation performed by the 
Remez exchange algorithm can be stated as follows: Find the 
set of filter coefficients that minimizes the maximum 
absolute value of the error, E(6), over the frequency range 
of interest. 

Eopt. = Min [Max 1E(8) 1] (2353) 
At this point, discussion of the weighting function, w(6), 
is in order. The purpose of this weighting function is to 
ensure a small tolerance for error in critical frequency 
ranges. 

If W(8) is large, this means a large deviation from 
the desired frequency response, Hp feu can not be toler- 
ated. Looking at equation (2.2), we see that if W(6@) is 
large, the difference between the desired and designed 
frequency response. [ Hp (e3®) - H (e38) ] has to be 
small to keep the weighted error small. Conversely, if W(8) 
is small, the difference [ Hp (e)®) - (eJ%) ] can be 
large and still meet the error criterion. Small values for 
W( 8) would be used in frequency bands where close approxi- 


mation to the desired frequency response is not critical. 
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2. Characteristics of the Remez Program 
The Remez exchange program requires the following 
Darameter specifications for a bandpass filter design. 
Frequencies are normalized with respect to the sampling 
Erecwenc y . For a bandpass filter the design procedure 
consists of specifying N the filter length, the left and 
meneeoee sGUEGrr f£Lrequency (Frp, Frp), the left and 


Mmeuclace stOport ~rreequency (Frc, Fro), and the ripple 


6 
ratio, K = oa which determines the desired weighting 
ie 
Ene elon wie> ) as: 
3 - ' 0S @ S 2HFr Os 
i ; 2TFRsSO8 TT 
Here 6, momene passband ripple, and 1s the stopband 


ripple, and AF is the transition band size. Figure 5 shows 
the frequency response of a bandpass filter and the 
frequency specifications. Since this optimal filter design 
procedure is based on the Chebyshev approximation which has 
ripples in the passband and stopband, the stopoff and cutoff 
Frequencies are dependent on the specified magnitude of the 
desired ripple. It is not possible to specify the 3-db 
frequency position for the designed filter directly in the 


Remez program. 


el 


1+, 


L\ /\I/\ 7\ 


ie 


ee ee en ae ae 
VS ee 


ETS One Frp Frs 


Figure 5 Frequency response of an optimal (minimax error) 
bandpass filter [Ref. 2]. 
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TO investigate the relationship between the filter 
length, N, and the transition band size, AF, the Remez 


program was used for a BPF with the following specifica- 


t HONS: 
Ec = O21 
Pap = 0.2 
Pep = 0.3 
ep = On4 
K ae 


Figures 6, 7, and 8, illustrate the frequency 
response obtained for N=31, 61, 81, respectively. The other 
specifications are kept the same. Looking at the left 
transition region we note, the longer Fondue fidter has a 


Sharper transition band. 


pO CO 
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Figure 6 The relation between filter length 
and transition band when N = 31. 
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Figure 7 The relation between filter length and 


transition band when N = 6l. 
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The relation between filter length and 


Figure 8 
transition band when N = 8l. 
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Figure 9 The relation between filter length and 
transition band slope. 


From the above study, it is known that if the 
filter length N is sufficiently large, the desired transi- 
tion slope can be achieved. The exact 3 db cutoff frequency 
of the designed passband is dependent on the filter length N 
and the stopband and passband specification. AS we can see, 
lt seems that there are relationships among the transition 
band (AF), filter length (N) and band rejection (B. R.) 
in FIR filter design using the Remez program. These will be 
discussed in Chapter III. 

Be STUB Y SEAM PLE : SPECIFICATIONS FOR A VIDEO INTERMEDIATE 
PREQUENGy (I. F.) FILTER 
The main effort of this section will be to take an 
analog filter and redesign the filter using the Remez 


IES) 


exchange technique. The first step is to convert the analog 
requirements to digital requirements using the formula 
ae QT, where T = 1/ sampling frequency and = analog 
frequency. There are advantages of using digital filters 
over analog filters. Namely, they are less likely to be 
affected by noise and are therefore more reliable than 
analog filters. Additionally, digital filters have many 
possible implementations, and the signal data can be stored 
without loss of accuracy. 
1. Analog Specifications 

The analog BPF commonly used for video signals has 
been selected as an example. The objective is to inves- 
tigate whether it is possible to implement a traditional 
analog BPF with digital circuits. In practice, two stages 
of I. F. amplifiers and the associated filters are used for 
processing of composite television signals. Figure 10 
lllustrates the transfer function curve of a full 4.25 MHz 
bandwidth amplifier. As can be seen it meets the following 


specifications. [Ref. 5]: 
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Figure 10 Typical I. F. amplitude response curve 
for full 4.25-MHz bandwidth. 
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a. Over-all Gain 

Tiemevera lime ain Gem tae yideo I-F. amplifier is 
determined based on the noise and Signal voltage from the 
meer Stage. The total gain in the picture I. F. amplifier 
muse be 12,500 times in voltage or 82 db. If we use a 
digital filter instead of the original analog amplifier and 
filters, this means that the band rejection ratio should be 
mee? , 500, thus making the rejection ratio less than 20 log 
ery r0, 000) = -80 db. 

b. Gain-bandwidth Product 

The gain-bandwidth product of a typical I.F. 
Stage 91S 53. This means that a gain of about 13 can be 
obtained at a nominal bandwidth of 4 MHz, or a gain of 17.5 
at 3 MHZ, or a gain of 17.5 at 3 MHz bandwidth. it sehws 
appears that four stages are necessary to achieve I. F. 
gains above 10,000 with an over-all bandwidth of 4 MHz. 
However, recent improvements in I.F. amplifier circuit 
arrangements have increased the gain-bandwidth product to 
about 75, making possible a gain of 10,000 in three stages 
at a full bandwidth of 4 MHz [{Ref. 7]. The present tendency 
is to design the I.F. amplifier for the maximum attainable 
bandwidth, consistent with the selectivity requirements, and 
thus assure maximum resolution in the image. In the design 
case considered here a passive digital filter is used to 
achieve the transfer function requirement. The gain- 


Sandwucdeneproduct 1S not.of concern. 


ey 


Sr Lee eas 

The third specification relates to the selec- 
tivity of the I. F: amplifler, @ind particularly scence 
crimination against three undesired carriers: the as- 
sociated sound carrier at 41.25MHZ whose total attenuation 
is -65 db, the lower adjacent channel sound I.F. carrier at 
47.25 MHz whose total attenuation is -20 db, and the upper 
adjacent channel picture I.F. carrier at 45.75 MHz whose 
total attenuation is from -10 to =-15 db. In other words, 


the analog I-F amplifier specifications are: 


1. Stop band 1: 0 (MHz) ------ 41.3 (MHz) 
Pass band 2: 41.73 (MHZ) ------- 45 (MHz) 
Stop panama: 46.5 (MHZ) -------- 100 (MHz) 

2. Band rejection ratio: t/12,500>= Geab, -80"d. 


3. Phase delay: It 1S required to be linear in the 
passband on both sides of the cutoff frequency from 
-0.87 to 0.87 radians. Use of an FIR filter will 
guarantee linear phase and constant delay for the 
frequency component. 
2. Translated Digital Specification 
To translate the analog Specifications ~ seme 
following are determined: 


1. Nyquist sampling frequency: F = 100 (MHz) 


This 1S an arbitrary choice dictated by forthcoming 


technology. 
2. Band 1: 6, = (2m x 0)/(100 x 106) = 0 
@2 = (2m x 41.3 x 106)/(100 x 106) = 0.8265 
Band 2: 63 = (2n x 41.73 x 106)/(100 x 106) = 0.834n 
64 = (20 x 45 x 10°)7 Goo 00) = oor 
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Baneeecee5 = (27 x—4ioe5 x 109 /(100 x 106) = 0,93, 


cee; Om uoeoy/ (100 x 105) = > 


For convenience, we normalized them to ™ radians. 


Band 1: 6,1 = 0 eS = of) 
maven T= y0mS a => = 6> = 0.413 

Bana 2: 0.834m/Tm = 03/0.5 => 93 = 0.417 
0.9m/T7 = ay 0 5 > 94 = 0.45 

Banda3 : O93 / = 0670.5 7 85 = 0.465 
1/1 = Ons => 66 =0.5 


3. Band Rejection Ratio and Phase Delay 

The band rejection is chosen based on the design. 
For a four stage implementation, each individual filter 
Stage rejection ratio should be at least -20 db. If the 
three stage method is chosen, then each individual filter 
rejection should be at least -26.64 db, (10,000 = 104 = 
21.52, 20 1Ogm2i.5 = 26.64dB). If two stages are chosen, 
any one rejection should be at least 40 db. 

As mentioned before, the phase delay requirement 1s 
automatically satisfied using FIR filters. Present Analog 
to Digital (A/D) converter technology can provide a maximum 
sampling frequency on the order of 50 MHz. The 100 MHZ 
sampling frequency desired iS not available at this time, 
however, it may soon be possible in the future. Although it 
is impossible to use this design in digital video processing 


now because of the high sampling rate requirement (100 MHz), 


ao 


the objectives of using this design illustration is to 


develop the FIR digital filter design methodology. 
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lita BANDPASS PILTER DESIGN 


The main purpose of this chapter is to use the analog 
Specifications of an I.F. filter as mentioned before in 
Chapter II and design a digital bandpass filter. To achieve 
this goal, several different methods are tried, i.e., single 
stage and multiple stage methods are considered for this 


design. 


A. SINGLE STAGE FILTER DESIGN 
1. Experimental design 
Prior to using the Remez program for filter design 
purposes, it is necessary to understand what band rejection 
characteristics are exhibited by the filter designed using 
this program. In other words, what effect does varying the 
program parameters have on the resultant filter characteris- 
Eee Ss 2 Lowmachievye —tmls = tamiliagization, the following 
experimental filters are used to show the characteristics 
of the design using the Remez program. 
Filter 1 is a given bandpass filter with transition 
Deiac = Oele, The normalized sampling frequency is 0.5. The 


specifications of this filter are as follows: 


stop band 1: 0.0 © Ore k8 
pass band: OR 8s EO Ors 2 
Stop band 2: O.42 5 GO Ou5 


ZA 


Filter 2 is another given bandpass filter with 


transition band=0.012 The speciftveat vensecre: 


Stopes band wie Oreo EO eee 
pass’ band: O22 cae Cras 
stop band 2: C357 tO Ors 


Filter 3 1s a bandpass filter with transition 


band=0.05. The specifications are as follows: 


stop “banda. se tO O23 
pass band: io to O23e2 
StOD ee band 2: Oo to Cm 


These three filters are then designed using the Remez 
program with different filter lengths, N, and weights. 
Table 3.1 contains a Summary of the results obtained from 
these experiments. Different N's are chosen because it is 
interesting to reveal the effect of varying the filter 
epoehie lot Two dr fterent sets Of ewergnts, Wy anda open 
chosen to control the passband and the stopband ripple. The 
transition band between stopband and passband (AF = cutoff 
frequency - stopoff frequency) is also varied. The purpose 
here 1s to examine the interdependency of the three 
parameters: AF, N and Band Rejection (B. R.). The relations 
are just like a triangle shown in Figure 11, but what are 
the exact dependencies among these three parameters? 

From Table 3.1, the limits of the band rejection with 
respect to the filter length can be obtained and are 


summarized in Table 3.2. 


Ze 


Figure 11 Filter length N, transition band AF, and 
achieved band rejection (B.R.). 


Pasice Se CHARACTERISTIC S@MOF FIMBEERS IN EXPERIMENT 
(SAMPLING FREQUENCY=100MHz) 


N (Filter W (Weighting) AF (Transition Band 
length) Rejection 





TABLE 3.2 CHARACTERISTICS OF FILTERS FOR A LIM ie 


REJECT LONM@r @20dr- 


W (Weighting) (Transition N (Filter Limit 
length) rejectio 
W> (pass (dB) 





From the investigation of how the program parameters 


affect filter characteristics, the following comments are in 


order: 


ie 


As shown in Table 3.1, when the W (Weighting function) 
and N (filter length) are fixed, the band rejection is 
proportional €o the AF (transition bandjmee(That i 
when AF becomes small, the rejection becomes small 
and poor). 


When the AF (transition band) and W (Weighting 
function) are fixed, the band rejection is propor- 
tional to the N (filter length). If more band 
rejection is desired, an increase in N, the length of 
the filter, is required. 


When W (weighting function) is fixed to satisfy a 
certain band rejection, the AF is inversely propor- 
tional to the N (filter length). From Table 3.1 we 
see that in order to satisfy the analog specifications 
with a band rejection of -80 db as discussed in the 
previous chapter, it is necessary to choose a filter 
length, N, ranging from 35 to 41 where the weighting 
function wy=w2=l. 
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ee 


Actual Bandpass Filter Design 


The digital specifications of the video I.F filter 


in Chapter II will be used, however, 


there 1s a finite limit 


on the filter length that we can use due to the characteris- 


tics of the Remez program and 


in mind these limitations, 


response illustrated in Figure 12 


hardware constraints. Bearing 
a filter with the frequency 


is obtained. Tis figkter 


has the following specifications: 


Peer length tue 11 


band specifications: 


Stop band 1: Ons LO) e413 

pass band: 0.417 to 0.45 

Stoo band 2: On46G5—8 to 0.5 
weightung Gaetorns: Il, 10, @ 


From Figure 12 the response has a dip of about -1 db in the 


passband and a band rejection 


the desired specifications 


requires a larger filter length. 


required filter length is: 


> 


ene | (4 7/N) 


of only -3 db. The reason why 


Cannot be achieved is that it 


According to Ref. 4, the 


[weceome | atamgs > are Stoport and cutoff frequency. 


Senee, Orda iney.415 < 0.465-0.45, 
peescclceumene smaller Fs - Fy, 
ante ey 0. 413 and 
0.004 2 (4n/N) 
N 21000 
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The filter length must be “greater St han 0 i a - ae 
impossible to implement using just one digital filter. 


Therefore it is necessary to pursue another design method. 


© 
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Figure 12 Single stage filter design of the 
TV I. Ff. Specifications: 
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DeeeULlI Phe STAGE FILTER DESIGN 
1. Cascade Implementation 

Using the associative and distributive properties, 
the equivalent unit sample responses for parallel and 
cascade implementation are shown in Figure 13. The question 
is which implementation method can achieve the desired 
specifications? CAD tools will be used to find our answer. 
The procedure is as follows. Using the Remez program the 
cascade response of the two filters is obtained. 

1. For example, we select a fourth order bandpass filter 
VeemmocmiamomweclfiCations: O., 0.1, 0.2, 0.3, 0.4, 0.5 
and filter 1's impulse response is obtained. 

eee oeommag another fourth order bandpass filter but 
shifting the band specifications to: Oreo, OS, 
awe, BO. 5Sye we Obtain filter 2's impulse 
response. 

3. Multiplying the two filters together in the frequency 
domain yields the overall frequency response of Figure 
14. 

The CAD tools used also allow us to enter the 
eorpeneanwtilter transfer Eunctions in polynomial form 
directly, and yield the response shown in Figure 15. After 
comparing Figure 14 and Figure 15, the procedure to cascade 
two filters 1s confirmed to be correct. We now proceed with 
I-F filter design using the Remez program. 

The issue now is to find the shift between the two 
component filters used in the cascade form implementation of 


Pie wrt icter wn. Figure 2. Assuming that a single 21 point 


filter is generated as the first filter and that the second 
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Figure 13 Equivalent systems in parallel and 
cascade forms. [Ref. 3] 
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Figure 14 Using FFT to process and get the combined 
result in cascade form. 
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mimeeter 1S the shifted version of the first one, we examine 


three cases. 


Peer. . 1] 
fay &f=0.05 
Poet =O O01 


ae * 
Te 
Ran, 


whe lal > See eer 
Pa 
0 nee, 


u 
3 


Figure 15 
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The three sets of band specifications are: 


CAO Uomo oe, O50, 0.40, 0.50 


OE Otero rere arue 50, 0.35, 0.50 


UPOp eo eto eo, U.o), 0.31, 0.50 


a) 2.0) a0 4.0 0) 
Frequency x10/7 Hz 


Using polynomial entry command to get 


frequency response in cascade form. 


Zo 


With a fixed weighting 


of 10, 1, and 10 several filters were 


designed and are summarized in Table 3.3. The frequency 


response curves of these same filters are illustrated in 


Figures 16 ‘tem3Gn 


SUMMARY OF CAS 


TABLE os. 3 


CADE FILTER DESIGN EXPERIMEN, 


Distance Band Extracted 
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length transition : 
2, ae a 


(normal freq.) Rejection from 
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Figure 16 is the original bandpass filter with a transition 
Penamorewl-lss The Specification for the frequency bands are: 
Ssuee Damd 1: 0.0 0.08 
pass band: O218 O...3:0 
stop band 2: 0.40 035.0 
Figure 17 is a shifted version of the original bandpass 
Peucemmoyemo.O2 of the normalized frequency. The specifica- 
tions for the frequency bands are: 
Steepeband 1: Q0.0 é 
pass band: OreZ OioZ 
stop band 2: 0.42 O50 
Figure 18 shows the response of a combined cascaded filter 
using the filters of Figure 16 and Figure 17. The observed 
bermicmerevectloms ratio in Figure 18 is about -80 db. The 
rejection ratios are listed in Table 3.3 along the row with 
Pee—-eeb,e ne =] O. 1, and shifted distance = 0.02. The 
rejection ratios corresponding to the cascade designs 


illustrated in the remaining figures are listed in the 


second to the last column of the table. 
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Figure 17 Shifted filter by 0.02 of the 
normalized frequency. 
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Figure 18 Frequency Response for the cascade realization 
Cree aoe Uctaatecmumertoure 16G0and Figure 17. 
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Figure 19 Shifted filter by 0.04 of the 
normalized frequency. 
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Figure 20 Cascade realization of Figure 16 and 
Figure 19 filters. 
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Figure 21 Shifted filter by 0.06 of the 
normalized frequency. 
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Figure 22 Cascade realization of the original filter of 
Pere MGwancdmenecr tl teme lor Figure 21. 
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Figure 25 Cascade realization of the original filter in 
Pigure 23 am@stme filter of Figure 24. 
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Figure 27 Cascade realization for the original filter in 
Figure 23 and the filter of Figure 26. 
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Figure 28 Shifted filter by 0.1 of the 
normalized frequency. 
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Figure 29 Cascade realization of original filter in 
Bidgirew2> andsemeuri1lter of Figure 28. 
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normalized frequency. 
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Figure 32 Cascade realization for the original filter 
in Figure 30 and the filter of Figure 31. 
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Figure 33 Shifted filtemis7 0c merc 
normalized frequency. 
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Figure 34 Cascade realization for the original 
filter in Figure 30 and the filter of Figure 33. 
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Figure 35 Shifted filter sy sles cee ene 
normalized frequency. 
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Figure 36 Cascade realization for the original 
filter in Figure 30 and the filter of Figure 35. 


From the preceding pages, where the cascaded 
Been -—“meaue chown In fF foures 18, 20, 22, 24, 26, 28, 30, 
32, 34, and 36, it appears that to obtain an increase in 
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bandwidth through the use of cascade filter realization is 
difficult, however, band rejection characteriseres mane 
improved using the cascade implementation. Thus we see that 
achieving the desired wide and flat passband response 
requires an alternative design approach. 

With analog active amplifier filters it is possible 
to stagger a filter with the shifted one, anda flat and 
Wide passband response can be obtained as in Figure 37. 
Figure 38 illustrates two FIR filters overlapped according 
to the procedure suggested in Figure 37. Figure wed 
illustrates their cascade combination, which is not wider 
than the original passive digital filter. The Remez program 
for passive digital filters 1s thus not able to obtain a 
wider frequency response, aS was the case with analog 
filters: 


ee a ge Da 
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Figure 37 Stagger-tuned analog design in 
cascade form. [Ref. 3] 
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Figure 38 Original filter and the shifted filter design. 
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Figure 39 Cascade filter combination. 


2. Parallel Implementation 
As mentioned at the beginning of this chapter, we 
can attempt to design our I-F filter using a parallel 
implementation. The following band specifications and a 


filter lengths of 21 are used. 


Stop band 1:0 .0 to Om 367 I 
pass band : 0.34071 to 0.3842047 
stop band 2: 0.4882047 to 035 


The resultant filter frequency response is illustrated in 
Figure 40. The filter is then shifted to the left so that a 
crossover occurs at magnitude 0.5 as shown in Figure 41. 


Using addition to combine the two filters together, the 


ie 


resultant frequency response is depicted in Figure 42. From 
this figure, it can be seen that the resultant response is 
wider, as we desired. 

Sencwmeaqhestionscmremaim,. such aS, 1S it possible to 
obtain the same frequency response using a direct implemen- 
Gatien Of a single filter transfer function? Which design 
method is best in terms of hardware implementation? These 


problems remain to be investigated. 
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Figure 40 Original filter for parallel implementation. 
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Figure 41 Shifted filter. 
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Figure 42 Combination filter in parallel implement. 
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iy. DECOMPOSITION tr oR IMPLEMENTATION 


This chapter will investigate the use of the zZero-one 
coefficient decomposition method to obtain any order filter. 
First, the brute force enumeration method is used to explore 
the possibility of taking an existing filter transfer 
function, and finding all possible combinations of decom- 
emosea factors. for the filter polynomial. The frequency 
responses are then obtained, and studied to find a relation- 
ship between the original transfer function and the 
decomposed transfer functions. Going through these results 
to get all the numerous errors lead us to an approach to 
find the best decomposed filter and save the number of 
multipliers. From these studies a computer flow chart is 
proposed where the best 0,1 decomposed filter design can be 


selected. 


A. HARDWARE REQUIREMENT 

As seen in Chapter III, uSing various multi-stage 
methods to design a filter with a wide flat passband is 
difficult. The Remez program has a limit in 1ts ability to 
produce filter transfer function patterns Similar to the 
Original analog filter design of Figure 37. It was also 
shown, using parallel implementation techniques, that good 


results like that shown in Figure 42 can be produced. 


> 


Now, our attention 1S drawn to implementation of the 
desired filter using direct design techniques and comparing 
the results with those obtained using parallel implementa- 
fon. Figure 42 and Table 3.1 were consulted to determine 
the direct design specifications: 0.07, 0-125, User 
0.3541, 0.4779, 0.5 and filter length = 41. The results are 
illustrated in Figure 43. Comparing these results with the 
filter of Figure 42, we note the direct design is better. 
However, a true comparison of the two techniques requires a 


comparison of the hardware implementation of the two 


designs. 
1. Direct Design 
41 = 
BO Z)uc= et a 
n=l 


From the above equation, the direct calculation of 
the FIR requires the following computations: 
Number of complex multiplications: 41 
Number of complex additions: 40 
Number of total delays: 40 
The hardware implementation is illustrated in Figure 44. 


2. For the Two Filters in Parallel Combination 


NO 
om 


Ze 
h(n) 2-n + § oh (nj) 271 
iL n=1 


H (2) = 


od 


We can get the following: 


complex multiplications = 42 


aC 


complex additions = 40 


delays = 20 


The hardware implementation is shown in Figure 45. From 


this comparison, 


the parallel implementation is shown to be 
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Figure 43 Direct Design Implementation. 





Figure 44 Direct implementation for a filter 
of length = 41. 
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Slightly advantageous in terms of total delay time, which is 
not a critical factor in most desmams. Thus, we can 
conclude that the direct implementation is superior to both 


cascade and parallel methods, because it requires less time 





for design and exhibits optimal Frequency response 


characteristics. 
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Figure 45 Parallel implementation for a filter 
of length = 21. 


B. ZERO-ONE DECOMPOSITION 

As we have seen direct filter design using the Remez 
program 1S a convenient tool. [Ref. 2] The question 1s, is 
there any method to decompose the desired filter transfer 
function using only 0, +1, -1 coefficients to save multi- 
plies? [Ref. 3] To investigate this question a brute force 
enumeration method will be used to approximate the desired 
filter. (There is no known methodology for decomposition.) 
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In this chapter a method is proposed to take advantage of 
tiemavdilable CAD tools to find a good decomposed filter 
transfer function that can approximate the desired filter. 
A fourth order transfer function was selected as our 
example: 
H(Z) = -0.08705 + 0.0128592-1 + 0.12619272 
+ 0.012859273 ~- 0.087105z74 [ae t) 
The frequency response is plotted in Figure 46 and a direct 
hardware implementation is shown in Figure 47 where 


multipliers are needed. 
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Figure 46 Frequency response for the transfer 
Punmcelom or: equation (4.1) . 
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The total possible number of factors are: 


factor l forme 9fi + aZ7) e%pge2 hee +a7eoe 
the possible number= 33 x 2 =54 


factor 2 form: (let admit © (ogee Cc 7mm encore 
the possible number= 32 x 22 =36 


(1 + aZ7l + bz-2) (1 + cz71 + qZ-2) , 
the possible number= 3¢ x 22 =36 


factor 3 Eormse. = a7 el + ia ana + Cy Lae dZ7*) 
the possible number= 3 x 24 =12 


factor 49form-( az 7 nee oie lees Clee + die 
the possible number= 24 =16 


Total number= 166. 


where a=b=c=d=+1,-1 or 0. 


0.012859 0.087105 





ee yt 


Figure 47 Direct hardware implementation 
of equation (4.1). 
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The total number of possible decompositions as enumerated 
1s 166, while the number of distinct decompositions is equal 
to 148. The difference is due to the fact that when factor 
forms 2, 3, and 4 are expanded, some are identical. The 
Bessie le cecomposition factors Of this example filter 
transfer function are: 

(L - 27l + 272) (2 + 27h + 772) = 1 + 272 + 274 

(L + 27-1 + z-2 + Z-3)(1 - z7t) = 1 - 274 

(1 - 271 + 277-2 - 27-3) (1 + 271) = 1 - 274 

(1 + 271 - z7-2 - 7-3) (1 - 27-1) = 1 - 2272 +274 

(1 + Jat) (1 + Tian) ie + ee ee a 

(1 + 27-1) (1 - 273) = 1 + g7l - 77-3 - 774 

(1 - z71)(1 + 273) = 1 - z7l + 77-3 - 7274 

fiecwoy to = Ze-) =-1 - 7 = 773 + 274 

(1 + 272) (1 - 272) = 1 - 274 


feline) 


it 
tr 
- 
WN 
i 
tr 
( 
N 
i 
WJ 
( 
N 
i 
> 


A Zet + 252) (1 - 274) 
(1 - z7~l + 272) (1 - 272) = 1 - 271 + 273 - 274 
ieee me Zee) (1 - gel + Z72) = 1 + 2-2 + 27% 
(L + 272 = 27-2) (1 + 272) = 1 + 271 + 273 + 274 
(1 - 271 - 272) (1 + 272) = 1 - 271 = 273 + 274 
(1 - 271 + 772) (1 + Zrl + 27-2) = 1 + 272 + 274 


(1 + 27h + 272 + 273)(1 + 27t) = 1 - 274 
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-10.00 aa 


Using CAD tools, the frequency responses of these decomposed 
approximation filters were plotted, and examined. Figure 48 
shows the best filter resembling the original filter whose 
decomposed transfer function is: 

H (Z) = 1 - 272 + 274 (4.2) 
Figure 49 shows the hardware implementation of this 
decomposed result. Comparing the decomposed expression in 
equation (422)—wieh themworl gana equation (4.1), it is 
clear that 2 less multiplies and 2 less adders are required. 
Thus, zero-one decomposition saves a lot of hardware in the 


implementation. 
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Figure 48 Decomposed filter frequency response that most 
closely resembles original filter frequency response. 
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Figure 49 Coefficient 0, +1, -1 decomposed hardware 
implementation of the filter in Figure 48. 


Besides the brute force enumeration method, is there any 
other method for obtaining the decomposed form of a desired 
filter? Is it possible to use CAD tools to find the answer? 
To answer these questions, a quantitative evaluation 
between Figure 46, the desired filter, and Figure 48, the 
designed filter, is necessary. As is well known, the 
minimum points in the frequency response plots of a filter 
indicate the places where the zeros of the filter transfer 
function occur, i.e., the roots of the numerator of the 
filter's transfer function. Thus, it's possible to plot the 
zeros of the original filter transfer function along with 
the zeros of the decomposed filter transfer function. This 


is shown in Figure 50. 
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It can be seen that there is a direct relationship 
between the distance of the zeros of the original and the 
decomposed filters and how well the decomposed filter's 
frequency response approximates that of the original. The 
relationship is that, the further apart the zeros are, the 
less accurate the approximating decomposed filter is. This 
measurement can be used as an error function to select the 
best decomposition approximation from several available 
ones for a particular filter transfer function. 

For example, the original transfer function equation 
(4.1) can be divided into four factored forms as follows: 

(Zt0289Z25 +) 0saol Zou 

(2405 8925 —=7 0R 4s eZee, 

(730596644 10 2orolE, 

(et 9664 —) 0237 Gone. 
The best decomposed approximation in Equation (4.2) is found 
to have the following factors: 

(Z2+0.886+30.5), (Z2+0.886-j0.5), 

(Z2-0.866+4970.49) , {2-05 866-70. 4c 


A plot of these two sets of roots is shown in Figure 50. 
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root 





PigugcmoOmeRoots plot of Ehe»original and the 
decomposed approximating filter. 


The error, that is, the distance between the decomposed 
AppLommaeing tilter roots and the original filter roots, is 
determined as follows: 
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dy = (0.25714 + (0.966)4 - 1 


do = 0/7 

error l= di ado. = Ue eae 

a3 = | \fioreeere > 0.517 - 11 + 
(0.892547)2 + (0.451289)4 - 1] 

dq = Ter ee 

SE Male |e a2 — da "- da. = 5572 


The total error for this decomposed approximation filter 1s: 


erres 


app ] tedaeo wan 


filte 


1 + error 2 = 0.318091. This calculation technique is 


rs. The frequency response plots are in Figures 


to-62. 


TABLE 4.1 APPROXIMATE ERRORS FOR THE DECOMPOSED 
APPROXIMATION FILTER 















Decomposed Fig # Error evaluation 








i 60) a 1.2909 
1.13178 


0.78507 
ee 0.31809 
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the possible decomposed approximation 
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Figure 51 Decomposed filter for the error 2.4938. 
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Figure 52 Decomposed filter for the error 2.3579. 
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Figure 53 Decomposed filter for the error 2.175. 
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Figure 54 Decomposed filter for the error 2.1496. 
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Figure 55 Decomposed filter for the error 1.86178. 
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Figure 56 Decomposed filter for the error 1.74975. 
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Figure 57 Decomposed filter for the error 1.46814. 
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Figure 58 Decomposed filter for the error 1.38558. 
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PacUEcmague Decomposed filter for the error 1.3677. 
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Figure 61 Decomposed filter for the error 1.13178. 
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Figure 62 Decomposed filter for the error 0.78507. 
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Why is it necessary to establish this error evaluation? 
With a valid error evaluation it is possible to determine 
the best decomposed approximation filter by finding the one 
with the smallest error. This technique lends itself to 
computer automation, and a flowchart illustrating a possible 


computer algorithm is proposed in Figure 63. 
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Figure 63 Algorithm for finding the best decomposed 
Apeaoximaeron filter of coefficients 0,+1,-1. 


es 


The use of the computer to determine the best decomposition 
formula is far superior than the original brute force 
enumeration method. 

From the above study, it is shown that it is possible to 
evaluate the decomposed approximation filter and automate 
the search procedures. The actual implementation of this 
procedure is not within the scope of this thesis study, 
however, it is proposed that another CAD tool could be 
devised to obtain the best decomposed approximation filter 


for any order filter. 
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V. CONCLUDING REMARKS 


This thesis has concentrated on procedures in FIR 
Seogitalbetiiter design using CAD tools. We saw that the 
filter length, N, has to be large enough so that the transi- 
tion band will sharp enough for it to meet design specifica- 
mronS . The relationships among AE ({(tvansition band), N 
(filter length) and B. R. (band rejection) in the Remez 
exchange algorithm were also analyzed clearly, and it was 
determined optimizing one variable has a trade-off effect on 
the others. Specifically, when W (weighting function) and N 
(filter length) are fixed, the band rejection is propor- 
meonal “to the Memteransition Band). When the AF 
(transition band) and W (weighting function) are fixed, the 
band rejection is proportional to N (filter length ). While 
when W (weighting function) is fixed, to satisfy a certain 
band rejection, the AF is inversely proportional to the N 
(filter length). These results can be used to guide us in 
Pueune Lllter design. 

In order to achieve the goal of designing an equivalent 
FIR digital filter that corresponds to an analog I. F. 
filter, direct and multiple stage filters were examined. 
There are two methods in multiple stage design: cascade and 
parallel. In the cascade design, it is not possible to make 
the digital filter passband flat and wide enough, but the 


iS 


band rejection is improved. In the parallel design, it 
seems we can more closely approximate the desired charac- 
tCERISt1 Gls - Finally, it was shown that direct design is 
Still better than any one of the multiple stage design 
methods when considering hardware requirements. 

Another reason why direct filter design is better than 
using the multiple stage design is that, although the Remez 
exchange algorithm is a very good tool for the design of FIR 
digital filters, its Chebyshev approximation characteristics 
yields ripples in the passband and the stopband. The exact 
3 db cutoff frequency and stopoff frequency are dependent on 
the filter length N and other specifications. Also, for 
the multiple stage design, there are difficulties in setting 
the cutoff and stopoflt frequency. 

Finally, the coefficient zero-~one decomposition 
approximation method was presented as a means to decompose 
any desired filter transfer function into a manner that 
saves multiplier hardware. In order to do that, the bimumee 
force enumeration method was used to investigate the 
problem. A simple filter was used as an example to search 
for the best approximation. Additionally an approach based 
on error evaluation between original filter and decomposed 
filter was proposed. In summary, it is possible to conceive 
an algorithm that will search for a good decomposed 
approximation filter, an idea that has potential for future 


work. 
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